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Abstract 

Faddeev calculations have been performed for nucleon-deuteron photodisin- 
tegration of ^He (^H) and proton-deuteron radiative capture. The bulk of the 
results is based on the AV18 nucleon-nucleon force and the Urbana IX three- 
nucleon force together with explicit exchange currents or applying the Siegert 
approach. Three- nucleon force effects are predicted for both processes and 
are qualitatively supported by available data. 
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I. INTRODUCTION 



Three-nucleon (3N) forces come more and more into the focus of few-nucleon studies. 
Pure 3N continuum measurements at the accelerator facihties lUCF [1], KVI [2], RIKEN [3] 
and RCNP [4] are performed around 100-200 MeV nucleon laboratory energies with the aim 
to confront data to theoretical predictions based on modern high-precision nucleon- nucleon 
(NN) forces only [5]. Clear-cut discrepancies for certain 3N observables against all those 
predictions can be considered to be good candidates for 3N force (3NF) effects. Thereby 
the theoretical investigations are based on numerically precise solutions of the 3N Faddeev 
equations. Then adding present day 3NF models and comparing to those data one tries 
to explore their strength and spin structure [6]. Right now these 3NF models are the 27r- 
exchange Tucson Melbourne (TM) [7], a modified version thereof, TM' [8], which is closer 
to chiral symmetry, and the Urbana IX [9] forces. Another path to learn about 3NF's is the 
study of the low lying spectra of light nuclei as performed in Greens function Monte Carlo 
calculations [10]. The inclusion of 37r-exchange ring diagrams with intermediate A's on top 
of the Urbana IX 3NF appears to be rather promising to improve the theoretical description 
of the spectra [11]. In all those investigations there is clear evidence found that present day 
NN forces alone fail to describe many of the studied observables and adding the presently 
available three- nucleon force models moves theory in the right direction. 

A recent approach towards nuclear dynamics is based on chiral perturbation theory, 
which is closely linked to QCD and develops nuclear forces in a systematic and controlled 
manner [12]. In that scheme which treats multi- pion exchanges explicitly and incorporates 
short range processes in the form of contact forces of increasing chiral dimensions also 3N 
forces are predicted and this consistently to NN forces. In [13] it has been demonstrated that 
like in the conventional approaches mentioned above 3N forces are unavoidable to predict 
binding energies of three- and four- nucleon nuclei as well as to remove discrepancies in 
certain 3N scattering observables. 

Electromagnetically induced reactions in the 3N system should also show effects of 3NF's. 
Since both 3N bound and scattering states enter into the nuclear matrix elements for photon 
induced processes and both types of states are affected by 3N forces it would be surprising 
if the various response functions for these reactions would be unaffected. In principle just 
by the continuity equation 3N forces lead also unavoidably to 3N currents. It is a quan- 
titative question based on current choices of nuclear force models to reveal signatures by 
switching on and off 3N forces. If certain observables are linked to binding energies and if 
all modern NN forces including the most recent ones based on chiral perturbation theory 
are unable to predict the experimental bound state energies but the inclusion of 3N forces 
is, we leave it to the reader to decide whether the changes in those observables are called 3N 
force effects or just binding effects. Apparently under these circumstances both are tightly 
bound together. It is only with oversimplified toy model NN forces which do not describe 
the rich NN data set that possibly experimental bound state energies can be achieved. Con- 
clusions based on those models should be taken with caution. The search for three- nucleon 
force effects in electromagnetically induced processes has been started before. For recent 
references see [14-17]. It is the aim of this paper to investigate the nucleon-deuteron (Nd) 
photodisintegration of ^He and as well as the time reversed proton-deuteron (pd) capture 
process using modern NN forces and various 3NF models. 
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The single nucleon current operator is supplemented by exchange currents either in the 
form of the Sicgcrt approximation or by explicitly including MEC's of the tt- and p-like 
nature. The treatment is carried through nonrelativistically, though presumably some of 
the data that we analyze, require at least relativistic corrections. 

Two-body photodisintegration of ^He (^H) has a long history. Barbour and Phillips [18] 
found that the incorporation of the interacting 3N continuum is crucial for the understanding 
of that process. They solved the 3N Faddeev equations, at that time of course based on sim- 
ple finite rank forces. This was taken up again more consistently by Gibson and Lehman [19] 
treating the 3N bound state and the final 3N continuum on equal footing. More recently 
the Bonn group [20,21] used more modern NN forces represented in finite rank form. They 
analyzed quite a few data and pointed out a correlation of a certain cross section peak height 
with the triton binding energy, an issue which we shall also address but now in the context of 
3NF effects. All the work mentioned relied on the Siegert approximation. The current was 
restricted to the dominant Ei multipole [20] or to the Ei and E2 multipoles [21]. In a recent 
paper [22] a benchmark was set on the total 3N photodisintegration cross section. There 
two quite different approaches, the Faddeev one and a hyperspherical harmonic expansion 
together with the Lorentz integral transform method were compared to each other using 
AV18 together with Urbana IX and reached a very good agreement. This documents the 
technical maturity of advanced present day approaches. 

Also for the pd capture process many experimental and theoretical studies have been 
performed in the past. We refer to [21,23,24] for references. Specifically we want to point to 
the theoretical investigations by Schadow et al. [21], Fonseca and Lehman [25] and to recent 
studies at very low energies by Viviani et al. [26]. 

The present investigation is restricted to nucleon- deuteron fragmentations in relation to 
^He (^H) and we refer to a forthcoming study for 3N fragmentations. 

The paper is organized as follows. In Sec. II we briefiy review our formalism and the 
dynamical ingredients. Our results for Nd photodisintegration are presented in Sec. Ill 
together with available data. The pd capture observables are discussed in Sec. IV. We 
summarize in Sec. V. 

II. FORMALISM 

We evaluated photodisintegration and pd capture before in [24] and [27,28], always using 

Faddeev-like integral equations. In [28] we formulated pd capture based on NN forces alone. 
There the Faddccv-like integral equation is identical to the one for Nd scattering. This is 
because in pd capture the 3N scattering state | ^i'^^^) enters directly. On the other hand 
in ^He photodisintegration the 3N scattering state (^(~) | is involved like in electrodisinte- 
gration of ^He. The way to derive Faddeev-hke integral equations in the latter cases is to 
apply the adjoint Moeller wave operator entering the nuclear matrix element to the right, 
namely onto the electromagnetic current operator and the ^He bound state [29]. This has 
the very big advantage that the driving term of that Faddeev-type integral equation is fully 
connected, namely proportional to the ^He bound state. Because of the formal identity 
of the nuclear matrix elements for photodisintegration and electron induced processes the 
same Faddeev-like integral equation is applicable. In the two cases only the components of 
the current operator in the driving term have to be chosen appropriately. Now using time 
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reversal symmetry one can relate the matrix elements for pd capture and photodisintegra- 
tion, which are evaluated quite differently. This is a highly nontrivial numerical test for the 
various complex numerical steps involved. 

In [24] we added a 3NF in the evaluation of the pd capture process and applied it to cross 
sections and several spin observables. The formalism was straightforward since we could use 
directly the Faddeev-like integral equation for Nd scattering including 3N forces as derived 
in [30]. In the same paper [24] we compared Siegert approximation to the explicit use of a 
restricted but possibly dominant set of mesonic exchange currents. 

Now for two- and three-body photodisintegration we would also like to formulate an 
extension including 3N forces. We shall proceed as follows. The nuclear matrix element for 
Nd-photodisintegration of a 3N bound state has the following form 

iV,^<i=(4-)|j.(g)|*bound), (1) 

where q is the asymptotic relative momentum between the proton and the deuteron, and 
jriQ) is the component of the 3N current operator. The scattering state (^^~'* | can be 
Faddeev decomposed 

(4"M=(4"M(1 + ^^), (2) 

where P, according to our standard notation [31], is the sum of a cyclical and anticyclical 
permutation. The Faddeev amphtude {tp'f^ \ obeys the Faddeev equation [30] 

(4"^ 1= ('^«- I I (^^^0 + (1 + P)vl'^Go{tGo + 1)) . (3) 

Here the channel state (^g- | enters, which is a product of a deuteron wave function and a 
momentum eigenstate of the remaining nucleon, the NN t-operator t acting on nucleons 2 
and 3, the free 3N propagator Go, the permutation operator P and v}^\ the part of a 3N 
force, which singles out particle 1. For our notation see [31] and [5]. 

Using Eqs. (1), (2) and (3) we can write the nuclear matrix element as 

N^'' = (0,- I (1 - K)-\l + P)j,{Q) I Abound), (4) 

where K is the kernel of the integral equation (3). We introduce 

\U) = {1- K)-\l + P)j,{Q) I Abound) (5) 

or explicitly the integral equation 

I [/) = (1 + P)jV(Q) I Abound) + {PtG^ + (1 + P)vt^G^{tG^ + 1)) I U). (6) 

This form is not yet suitable for numerical applications because of the presence of P to the 
very left. This has already been noted at the very beginning of our numerical 3N studies 
using nuclear forces without finite rank representations [32]. To rewrite Eq. (6) into a 
suitable form we use the following obvious identities 

(1 + P) = ip(l + P) (7) 
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-P(P -1) = 1 (8) 

and obtain 

{P-1)\U) = {P- 1)(1 + P)MQ) I Abound) + 

(P - 1)P (tGo + ^(1 + P)vi'^Go{tGo + 1)) ip(P - 1) I t/). (9) 

This Faddeev-like integral equation is suitable for numerical implementations and has the 
form 

I ^) = (1 + P)jV(Q) I Abound) + 

(tGoP + ^(1 + P)v}'''Go{tGo + 1)P) I U) (10) 

with 

\U)^(P-1)\U). (11) 
Then the nuclear matrix element results as 

N^'-l{<P,^\P\U). (12) 

In view of a forthcoming paper we also describe now the treatment of the complete 3N 
break-up process. The nuclear matrix element is 

iVf=(4-i I J,(Q) Inbound). (13) 

The asymptotic momenta of the three nucleons are given by standard Jacobi momenta p 
and g [31]. The Faddeev amplitude corresponding to the scattering state in Eq. (13) is now 
defined via 

(4^1= ^-\pq \ (14) 

where K is the same kernel as used before in Eq. (3) and 

(-)(pg|= {<f>o\{tGo + l). (15) 

It is to be noted that the free two-body subsystem state in {(pQ \ is properly antisymmetrized. 
Here (0o | is the free 3N state. Following the same steps as above one ends up with 

ATf = i(0o|(tGo + l)P|C/), (16) 

where | U) is as given above. 

In the actual numerical calculation we used, however, another form, which for the purpose 
of completeness, we would also hke to present here. The reason for that is that at the time 
of the installation, the very heavy numerical tasks were more easily performed with already 
existing building blocks. That alternative forms are for the Nd-break up 
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N^" = (<^,- I (1 + P) I MQ) I Abound) + (<^,- \P\U') (17) 

and for the 3N break up 

AT^N = (0, I (1 + P)jV(Q) I Abound) + (00 I ^^0(1 + P)jr{Q) \ ^bound) + 

{(l)o\P\U') + {<Po\tGoP\U'). (18) 
The Faddeev-hke integral equation for | U') reads then 

I U') = (^tGo + ^(1 + P)Vl''^Go{tGo + 1)) (1 + P)jr{Q) I Abound) + 

l^tGoP + ^(1 + P)l^i'^G'o(tG'o + 1)P) I t/')- (19) 

The equivalence between the matrix elements (12) and (16) on the one hand and (17) 
and (18) on the other hand is demonstrated as follows. From Eqs. (10) and (19) we have 

I u)- I U') = 

(1 - K)-' ((1 + P)jr{Q) I Abound) - ^^0(1 + P)jr{Q) \ *bound)- 
^(1 + P)V}'^Go{tGo + 1)(1 + P)jr{Q) I Abound)) = 
(1 - K)-l (l - tGo - ^(1 + P)Vl'^Go{tGo + 1)) (1 + P)jr{Q) I Abound). (20) 

Using again Eq. (7) and the form of the kernel K this is 

I u)- I U') = 

(1 - K)-^^{P - 1 + 1 - X)(l + P)j,(g) I Abound) = 

^(1 - K)-\P - 1)(1 + P)UQ) I Abound) + ^(1 + P)jr{Q) I Abound) = 
^(1 - K)-\l + P)j,{Q) I Abound) + ^(1 + P)jr{Q) I Abound) = 

^|C/) + ^(l + P)j,(Q)|*bound). (21) 

Thus 

^ I ^)- I U') = i(l + P)jriQ) I Vl/bound). (22) 

Now it is a simple task to verify that the two expressions for the 3N break-up amplitude 

(16) and (18) arc identical. This is also true for the two Nd break-up amplitudes (12) and 

(17) . It is efficient to evaluate also the pd capture process using time reversal in terms of 
the formalism just described. This is what we do in this paper. 

III. RESULTS FOR PD (ND) PHOTODISINTEGRATION OF ^^HE (^H). 

We use the AV18 NN force [33] combined with the Urbana IX 3N force [9]. By construc- 
tion that force combination describes the binding energy correctly. It over binds however 
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slightly the ^He bound state energy by 21 keV [34]. For the convenience of the reader we cite 
the theoretical binding energies : -7.628 MeV for '^H and AV18 alone, -8.48 by construction 
including Urbana IX, -6.917 MeV for ■^He and AV18 together with the Coulomb force and 
-7.739 MeV including in addition Urbana IX. The latter value is to be compared to the ex- 
perimental value -7.718 MeV which is shghtly different, but this should be of no significance 
for the present studies. Unfortunately we are still unable to include the Coulomb force into 
the pd continuum, which causes an inconsistency of unknown magnitude. Of course at least 
for the higher energies studied in this paper we expect minor Coulomb force effects as it is 
supported by pure pd scattering investigations [5]. Our calculations are fully converged by 
choosing total two- nucleon angular momenta up to jmax=3 and total 3N angular momenta 
up to Jmax— 15/2. It turned out to be sufficient to keep the 3N force different from zero for 
total 3N angular momenta up to J= 7/2. The standard nonrelativistic form of the single 
nucleon current operator [35] is supplemented by exchange currents via the Siegert approx- 
imation. We use it in the form as detailed in [24]. In our treatment electric and magnetic 
multipoles are kept to a very high order (6-7) and no long wave length approximation is 
used. Both ingredients are important as has been shown for instance in [22]. The formal- 
ism is performed throughout in momentum space. As is well known the Siegert approach 
corrects for many body currents only in the electric multipoles. Available models for two- 
body currents should then be added for the magnetic multipoles. This, however, is not yet 
included in this work. On the other hand we also use explicit exchange currents of the tt- and 
p- like type consistent to the AVI 8 NN force. Again this is not yet a complete approach since 
further pieces in the AV18 NN force have no counterparts in two-body currents which would 
be also required to fulfill the continuity equation. This needs further investigations [36], 
though the expectations are that with the n- and p- like parts the dominant currents are 
taken into account. If the continuity equation would be fulfilled in relation to all parts of 
AVIS the Siegert approach with respect to the electric multipoles would be essentially equiv- 
alent to these explicit MEC's (except of additional 3N force effects included in the Siegert 
approach and less important terms of higher multipolarities, see [29]). Our aim here is not 
to forward the theory of the electromagnetic current operator but to apply what is often 
called " the standard model of nuclear physics" to the complex 2-body photodisintegration 
or pd capture processes, which has not been done before to the best of our knowledge. 

We show in Fig.l the angular distribution for pd photodisintegration of ^He against the 
angle between the outgoing deuteron and the incoming photon direction in the laboratory 
system. The photon energies vary between 10 and 140 MeV. At the two lower energies 
the cross section maximum is decreased by adding the 3NF. A related effect has been seen 
before in [20] using different NN forces. Thereby it was found, that with increasing binding 
energy the value of the maximum decreased. This can be considered as a scaling behavior 
with the 3N binding energy. It ceases to be valid for the higher energies, where the results 
including the 3N force overtake the ones without. At about E^= 28 MeV the 3N force 
effects for the process ^He(7, d)p in that observable vanish. In relation to that scaling at 
low energies one can ask the question whether the 3N force contributions in the continuum 
are critical for that result. To that aim we performed calculations where we switched off 
the 3N force in the continuum but kept it in the '^He bound state. Thereby the ^He binding 
energy did not change. (Note this is not a consistent calculation and necessarily induces 
spurious effects). The effect is a decrease of the cross section of about 16% at E^=100 MeV 
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in relation to the difference of the results when the 3N force is dropped totally. At =10 
McV the presence or absence of the 3NF in the continuum had only a very tiny effect. 
Nevertheless we would like to repeat that the correct binding energy could only be achieved 
by adding the 3NF to the current most modern NN forces. In [37] energy weighted sum 
rules for the A=3 photodisintegration cross sections based on the electric dipole operator 
have been investigated. They link the energy dependence of the cross section through the 
integrals to expectation values of the ground state wavefunctions, which are affected by the 
3N forces and consequently depend on the binding energy. It appears worthwhile to check 
the assumptions and approximations in [37] from the point of view of present day forces, 
wavefunctions and currents. 

Have these effects already been seen in some data? We are aware of cross section data at 
the deuteron laboratory and cm. angles of 90 ° as a function of [39,38]. They are shown 
in Fig. 2 in comparison to our theoretical results. We see the crossing of the theoretical 
curves without and with 3NF's around 25 MeV and indeed the data support the decrease 
of the cross section at lower energy values as predicted by including the 3NF. At higher 
energies the effects of the 3NF appear to be somewhat too strong in case of the Ticcioni 
et al. data [38]. It is possible that the overshooting of the theory at the lowest energies is 
partially due to the neglection of the pd Coulomb force effects in the continuum. Precise 
new data would be very welcome. 

For E^— 120 MeV and higher photon energies we are aware of another set of data for 
that process [40,41]. The deuteron laboratory angle is 103 ° now. In this case the addition 
of the 3NF is clearly supported by the data. This is shown in Fig. 3. In this case we used 
the explicit tt- and p- like MFC's since the energies are higher and Siegert as a low energy 
approximation is less suited. 

There are also total nd and pd cross section data for the processes ^H(7, and 
^He(7, d)p. We show them in Figs. 4 and 5 based on the Siegert approach. Clearly the 
old data for the nd cross section have too big error bars to be conclusive. In the case of the 
pd cross section the inclusion of the 3NF's deteriorates the agreement somewhat for the low 
photon energies. The nd cross section data have been displayed before in [22] . 

Summarizing, the comparison with the angular distribution data appears to be in qual- 
itative agreement. New improved data would be welcome and a more refined treatment of 
two-body currents is required. 

In order to provide information on the dependence of the cross sections on the choice of 
forces and currents we display in Table I results for the total two-body photodisintegration 
cross section of ^He (^H) at three energies. At 12 MeV we see 5 % (10 %) spreads with 
(without) 3NF's. At the higher energies the spreads are negligible which points to a certain 
stability of the results and helps to identify 3N force effects. Precise data, however, would 
be required. 

IV. PD CAPTURE CROSS SECTIONS 

In Ref. [24] cross sections and spin observables for pd capture have been investigated 
at proton laboratory energies Ep between 5 and 200 MeV (corresponding to deuteron labo- 
ratory energies E^ between 10 and 400 MeV). The emphasis was on testing the sensitivity 
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of pd capture observables to changes in the choice of NN forces and to compare the pre- 
dictions of the Siegert approximation to the ones including exphcitly vr- and p-hkc MEC's. 
Please note as pointed out above that both approaches in the way we treat them are ap- 
proximate and therefore the comparison is more of qualitative nature. We found that at 
low energies Siegert and MEC predictions are rather close together, whereas at the higher 
energies differences showed up. In the context of the Siegert approach the predictions based 
on different NN forces turned out to be rather close together, which is a satisfactory result, 
since it demonstrates stability. The agreement with the data was mostly good, but also 
clear discrepancies were present, which call for an improvement of the dynamical input. It 
is the aim of this paper to include 3NF's, which in the previous work were only marginally 
investigated. 

In the following we show our results for cross sections between Ed= 19.8 MeV and Ed= 
400 MeV. In all of the following figures four theoretical curves are displayed. They are 
based on the Siegert approach, the single nucleon current together with explicit tt- and p- 
hke MEC's consistent to the NN force, the Siegert with 3NF and MEC's with 3NF. Let us 
denote these four choices by ' and b', for short. 

We see in Fig. 6 the pd capture cross section at Ed = 19.8 MeV. In both cases, Siegert and 
explicit MEC's, the inclusion of the 3NF decreases the cross section; in case of Siegert the 
decrease is much stronger. The choices a', b and b' are well within the error bars and only a 
is significantly too high. At E^ — 95 MeV the cross section data are fairly well described by 
all four choices. This is displayed in Fig. 7. As already seen in the Nd photodisintegration 
the theoretical cross section increases by including 3NF's. This is in agreement with our 
findings for pd capture [24]. 

Finally we show the cross sections for Ep= 100, 150 and 200 MeV (corresponding to 
Ed^ 200, 300 and 400 MeV) in Fig. 8. The cases with the exphcit MEC's and 3N forces 
(b') describe the data best (except for small angles). The choice a clearly underpredicts the 
maxima. 

V. SUMMARY AND CONCLUSIONS 

We presented the formalisms for including 3NF's into the Faddeev framework for photo- 
disintegration of three nucleon bound states. The resulting equations are solved rigorously 
using high precision nuclear forces: AV18 together with the Urbana IX 3NF or CD Bonn 
with TM 3NF. Many-body currents are included either in the form of vr- and p-like ex- 
changes related to AV18 or via the Siegert approach where the latter corrects only electric 
multipoles for many body currents and the former does not include all two-body currents 
consistent to AV18 (in the sense of fulfiUing the continuity equation). Thus both ways of 
going beyond the single nucleon current are approximate but currently used in the literature. 
The calculations are nonrelativistic but employ state of the art dynamics. We posed several 
questions. How well do the Siegert approach and the explicit use of the tt- and p-like MEC's 
compare with each other ? Our results displayed for pd capture show differences between 
the two approaches which calls for improvements either by adding two-body currents for the 
magnetic multipoles in the Siegert approach or by adding at least the currents beyond the 
TT- and p-like parts required for the consistency to AV18. Qualitatively, however, the two 
approaches give similar results. 
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Another even more central question in that paper was to shed hght on possible 3NF 
effects. In case of pd photodisintcgration of ^He we compared theoretical predictions without 
and with a 3NF. We found a clear signature in adding the 3NF. The maxima are decreased 
at low energies and increased at high energies. The turning point is around E^=28 MeV. 
At the low energies this can be considered as a scahng effect with the ^He binding energy 
but one has to note that based on present day NN forces the '^Hc binding energy can only 
be achieved if a 3NF is added. These 3N force effects up to about 60 McV are too small to 
be verified by the presently available data. However at the lower energies about 10 MeV 
our theoretical predictions including 3N forces are clearly too high which might be due to 
the neglected Coulomb force in the continuum. At the higher energies > 120MeV the 
effects arc larger and qualitatively supported by the data. One should, however, be aware 
that beyond the 7r-threshold we certainly leave the validity of our nonrelativistic framework. 
In case of the pd capture at the higher energies, E^=100 MeV and above, explicit use of 
MEC's together with the 3N force model shows a tendency to move theory better towards 
the data than without 3N forces. The failure at the smaller angles shows however that 
some ingredients are missing. Overall we demonstrated that 3NF's can be incorporated into 
such a complex 3N reaction process and effects are visible related to the models used. An 
improved theoretical treatment of many-body currents and more precise data are needed to 
achieve a clearer view towards 3N force effects. 

Altogether the shifts caused by the Urbana IX 3NF on top of the AV18 NN force and 
explicit MEC's is supported by most of the data we analyzed. The Siegert approach is 
less successful. The use of other force combinations as exemplified in the total two-body 
photodisintcgration cross section does not lead to alarming variations. High quality data 
would be very helpful to challenge theory more strongly. 
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TABLES 





= 12 fMeVl 


= 40 [MeV] 


= 120 [MeV] 

/ L J 


AV18-Siegert 

AV18-MEC 

CD Bonn2000-Siegert 


1.086 (1.056) 
0.953 (0.949) 
0.997 (0.980) 


0.160 (0.168) 
0.156 (0.155) 
0.163 (0.169) 


0.016 (0.015) 
0.017 (0.015) 
0.017 (0.016) 


AV18+UrbanaIX-Siegert 

AV18+UrbanaIX-MEC 

CDBonn2000+TM'-Sicgert 


0.932 (0.882) 
0.934 (0.915) 
0.917 (0.889) 


0.173 (0.180) 
0.172 (0.169) 
0.170 (0.176) 


0.020 (0.018) 
0.020 (0.017) 
0.020 (0.018) 



TABLE I. The total cross section (in mb) for two-body photodisintegration of ^He (^H). 
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FIGURES 
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FIG. 1. Deuteron laboratory angular distribution for the process '^He(7, d)-p at different photon 
energies E^. Curves show results of calculations with the AV18 NN and Urbana IX 3NF forces 
(solid) and with the AV18 NN force alone (dashed). The current is treated in the Siegert approach. 
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FIG. 2. Deuteron angular distribution for the process '^He(7, d)p at the given deuteron angle 
as a function of the photon energy E^. Curves as in Fig. 1. Since the kinematical shift from the 
laboratory to the cm. system is not significant, we combine the data for the 90 ° laboratory angle 
(full dots with horizontal and vertical error bars [39]) with the ones for the 90 ° cm. angle (open 
circles [38]). 
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FIG. 3. Deuteron angular distribution for the process ^He(7, d)p at given laboratory angle as 
a function of the photon energy Ej. Curves show results of calculations with the AV18 NN and 
Urbana IX 3NF forces (solid) and with the AV18 NN force alone (dashed). Explicit tt- and p-like 
MEC's are included in the current operator. Data are from [40] (x-es) and [41] (circles). 
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FIG. 4. Total cross section for the process ^H(7, d)n as a function of the photon laboratory 
energy E^. Curves as in Fig. 1. Data are from [42] (x-es), [43] (squares), [44] (triangles). 
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FIG. 5. Total cross sections for the '^He(7, d)p process as a function of the photon laboratory 
energy E^. Curves as in Fig. 1. Data are from [45]. 
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FIG. 6. The photon angular distribution for pd capture at Efi= 19.8 MeV against the cm. 7-d 
scattering angle. The curves describe the Siegert (dashed-dotted) , the single nucleon plus MEC 
(dashed), Siegert with 3NF (dotted) and the single nucleon plus MEC with 3NF (solid). These 
four cases are called a, b, a' and b' in the text. Data are from [46]. 
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FIG. 7. The same as in Fig. 6 for Ed= 95 MeV. Data are from [47]. 
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FIG. 8. The photon angular distribution for pd capture at three different proton laboratory 
energies {Ep = 100, 150 and 200 MeV from above to below) against the cm. 7-p scattering angle. 
Curves as in Fig. 6. Data are from [48]. 
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